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Abstract. We have analyzed ACS/WFC and WFPC2 images from HST, as well as ground-based data to study the radial 
distribution of the double sub-giant branch (SGB) recently discovered in the Galactic globular cluster NGC 1851. We found 
that the SGB split can be followed all the way from the cluster center out to at least 8' from the center Beyond this distance out 
to the tidal radius at ~ 1 1.7 arcmin, there are simply too few SGB stars to identify the sequences. The number ratio of the bright 
SGB to the faint SGB stars shows no significant radial trend. Furthermore, we have found that the ratio of blue horizontal-branch 
(HB) stars plus RR Lyrae to the red HB stars also remains constant from the cluster center to the outer envelope. 



1. Introduction 

Precise Hubble Space Telescope (HST) photometry has pro- 
vided evidence that NGC 1851 hosts two distinct sub- 
populations of stars (Milone et al. 2008) as indicated by a clear 
bifurcation of the sub-giant branch (SGB) in its color magni- 
tude diagram (CMD). This discovery has sparked new interest 
in this object and, consequently, numerous efforts aimed at bet- 
ter understanding how this cluster formed and evolved. 

Milone et al. (2008) suggest that two star-formation 
episodes delayed by about one Gyr could explain the observed 
split of the SGB. As an alternative scenario, Cassisi et al. 
(2008) suggest that the SGB split can be explained by the pres- 
ence of two stellar populations, one with normal a-element en- 
hancement, and the other characterized by a peculiar CNONa 
chemical pattern with C-i-N+O abundance increased by a fac- 
tor of ~ 2. Two such populations could account for the observed 
bright and faint SGB (hereafter bSGB and fSGB), respectively, 
without requiring any significant age difference. 

Interestingly, this latter scenario seems to be supported by 
early spectroscopic measurements (Hesser et al. 1982) which 
indicate the presence of two groups: CN-strong and CN-weak 
stars. In addition, the work of Calamida et al. (2007) shows 
that in the Stromgren (mi, u - y) CMD, the red giant branch 
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(RGB) of NGC 1851 splits into two sequences that, Uke the 
SGB split, can be explained by two populations with different 
CN abundances. 

A more recent spectroscopic investigation by Yong & 
Grundahl (2008) determined the chemical composition for 
eight bright giants in NGC 1851. Their analysis revealed 
large star-to-star light-element-abundance differences of the el- 
ements Zr and La. These i-process elements are coiTelated with 
Al and anticorrelated with O. Furthermore, the Zr and La abun- 
dances appear to peak around two distinct values. Yong et al. 
(2009) show that the C-i-N-i-O abundance exhibits a large spread 
(~ 0.6 dex), giving further support to the Cassisi et al. (2008) 
scenario. They also found a correlation of Na, Al, Zr, La abun- 
dance with C-hNh-O, as expected in the scenario in which in- 
termediate mass AGB stars are responsible for globular-cluster 
light-element abundance variation. 

Another important peculiarity of NGC 1851 is that it is one 
of only a few examples of a bimodal horizontal branch (HB) 
cluster. On the basis of detailed numerical simulations, Salaris 
et al. (2008) found that it is possible to account simultaneously 
for the various empirical constraints such as the HB morphol- 
ogy, the star counts along the HB as well as the ratio between 
faint SGB stars and bright SGB ones, if all the upper SGB stars 
and a small fraction of lower SGB stars evolves into the red 
HB, while most of the lower SGB stars populates the blue HB 
(including the RR Lyrae variables). 

Therefore it is tempting to associate the bSGB stars with the 
CN-normal, s-process-normal stars and with the red HB, while 
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the fainter SGB should be populated by CN-strong, .s-process- 
element-enhanced stars which should evolve mainly onto the 
blue HB. This scenario implies that fSGB stars correspond to 
the second generation and that they formed from material pro- 
cessed through a first generation of stars. 

The study of the spatial distribution of the two stellar popu- 
lations associated with the double SGB occupies a pivotal role 
in the current research on multiple stellar populations in GCs. 
Milone et al. (2008) studied a CMD containing stars with radial 
distances smaller than ~2.5 arcmin from the cluster center. In 
this paper, we extend the analysis to the outermost regions of 
the cluster Our radial extension covers the whole cluster from 
the center to beyond the tidal radius (r,iij„i-llJ arcmin, Harris 
1996, 2003). 

A first attempt at constraining the radial distribution of the 
two stellar populations associated with the double SGB comes 
from the recent work by Zoccali et al. (2009). They analyzed 
VLT-FORS V, I images, of the South West quadrant of the clus- 
ter with the aim of following the extent of the double SGB from 
~1 to ~13 arcmin from the cluster center. Zoccali et al. (2009) 
claimed that the percentage of fSGB stars, which is ~45% in 
the innermost region decreases at ~ 1.5 arcmin from the center 
and completely disappears at ~ 2.4 arcmin. The more extensive 
and higher-precision data set presented below does not confirm 
the results of Zoccali et al. (2009). 

The plan of the paper is as follows. In the next section we 
describe the data and the data reduction. The CMDs from dif- 
ferent data sets are presented in Sect. [3] In Sect. |4] we describe 
the criteria used to define the NGC 1 85 1 sub-populations and 
how we measured their relative frequency at several radial dis- 
tances. We present the radial distribution of the two stellar pop- 
ulations in NGC 1851 and discuss our results in Sect.|5] 



Table 1. Description of the HST data sets used in this paper. 



INSTR 


DATE 


NxEXPTlME 


FILT 


PROGRAM (PI) 


ACS/WFC 


May 01 2006 


20s+5x350s 


F606W 


10775 ( Sarajedini) 


ACS/WFC 


May 01 2006 


20s+5x350s 


F814W 


10775 ( Saiajediiii) 


WFPC2 


Nov 07 2007 


I0x230s 


F450W 


lI233 (Piotto) 


WFPC2 


Nov 07 2007 


5x200s+5x230s 


F814W 


11233 (Piotto) 



2. Observation and data reduction 

In order to study the radial distribution of SGB and HB stars in 
NGC 1851, we considered three distinct data sets. 

To probe the most crowded regions of the cluster we took 
advantage of the high resolving power of HST, using im- 
ages collected with the Wide Field Channel (WFC) of the 
Advanced Camera for Survey (ACS) and with the Wide Field 
Planetary Camera 2 (WFPC2). NGC 1851 is not in a particu- 
larly dense region of the Galaxy (fl = 245°, b^^ = -35°), and 
ground-based observations can provide photometry outside the 
crowded central region with precision comparable to that of the 
more central regions observed with the HST cameras. 
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R.A. 

Fig. 1. For the innermost 8 arcmin, we show the footprints of 
the HST-images used in this paper superposed on a DSS image. 
In cyan ACS/WFC (circles highlight the corners of the first chip 
of ACS), in red WFPC2. Large circles show the different radial 
bins described in Section|4] The different colors refer to differ- 
ent data-sets: cyan for ACS/WFC, red for WFPC2, and green 
for ground-based data-set. 



A brief description of the HST images used in this work is 
given in Table [1] while Fig. [T] shows a finding-chart of their 
footprints. 

2.1. The ACS/WFC data-set 

The HST ACS/WFC images come from GO- 10775 (PI: 
Sarajedini, see also Sarajedini et al. 2007) and were presented 
in Milone et al. (2008); we have used the output of the reduc- 
tion described in Anderson et al. (2008). In brief, the proce- 
dure analyzes all the exposures of each cluster simultaneously 
to generate a single list of stars for each field. Stars are mea- 
sured independently in each image by using the best available 
PSF models from Anderson & King (2006). 

This routine was designed to work well in both crowded 
and uncrowded fields, and it is able to detect almost every star 
that can be perceived by eye. It takes advantage of the many in- 
dependent dithered pointings of each scene and the knowledge 
of the PSF to avoid including artifacts in the list. Calibration of 
ACS photometry into the Vega-mag system was performed fol- 
lowing recipes in Bedin et al. (2005) and using the zero points 
given in Sirianni et al. (2005). 

2.2. The WFPC2 data-set 

The HST WFPC2 images come from GO-11233 (PI: Piotto), 
a proposal specifically dedicated to photometric detection of 
multiple populations. 
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The WFPC2 images have been reduced following the 
method of Anderson & King (2000), which is based on 
effective-point-spread-function fitting. We corrected for the 
34"' row error in WFPC2 CCDs (see Anderson & King 1999 
for details) and used the best distortion solution available, as 
given by Anderson & King (2003). Photometric calibration has 
been done according to the Holtzman et al. (1995) Vega-mag 
flight system for WFPC2 camera. 

2.3. The ground-based data-set 

The ground-based data are taken from the image archive main- 
tained by one of us (Stetson 2000). The observations used 
here include 545 images from 14 observing runs with the Max 
Planck 2.2m telescope, the CTIO 4m, 1.5m, and 0.9m tele- 
scopes, and the Dutch 0.9m telescope on La Silla. Any given 
star may have as many as 69 independent measurements in the 
B filter, 78 in V, and 56 in /; among these, 62, 70, and 56, 
respectively, were taken on occasions that were judged to be 
of photometric quality. These data were reduced following the 
protocol outlined in some detail in Stetson (2005). We have 
complete photometric coverage of the cluster field out to a ra- 
dius of 14.5 arcmin, and partial coverage to 25 .9 arcmin, but we 
will restrict our discussion here to the area within 12.5 arcmin, 
which is slightly larger than the tidal radius of 11. 7 arcmin. 

3. The CMDs 

Figure |2] shows the CMDs from HST observations that cover 
the densest regions of the cluster. The left panel shows the 
CMD already published by Milone et al. (2008), and the right 
panel shows the newly derived CMD from the WFPC2 data 
described in Sect. 2.2. 

The CMDs from ground-based photometry are shown in 
Fig. |3] where we have rejected all the stars within 2.5 arcmin 
from the cluster center. The split of the SGB is clearly visible 
both in the (V - I) vs. V and the (B - I) vs. V CMD, and the 
split region is highlighted in the inset where we show a zoom 
of the SGB region. Note that the fact that we see the SGB split 
beyond ~ 2.5 arcminutes is already in disagreement with the 
claims by Zoccali et al. (2009) who analysed the S-W quadrant 
and found no evidence for the presence of fSGB stars at ra- 
dial distances larger than ~ 2.4 arcmin. The left panel of Fig. |4] 
shows our y vs (Z? - /) CMD from ground based data for stars 
with radial distance from the cluster center larger than 2.5 ar- 
cmin. We marked in red the stars that more likely belong to 
the fSGB. In the four right panels we plot the same CMD for 
stars in the four quadrants. The quoted numbers are the number 
of selected fSGB stars with Poisson errors. Because of the un- 
certainties we conclude that there is not significant difference in 
the distribution in the CMD of fSGB stars in the four quadrants. 
We detected 23 probable fSGB stars in the S-W field covered 
by Zoccali et al. (2009) most of them belonging to the lower 
part of the fSGB. Possibly, the difference between our result 
and the one obtained by Zoccali et al (2009) comes from the 
small number of fSGB stars with radial distance greater than 
~2.5 arcmin in the S W quadrant. It must be noted that the CMD 
of Zoccali et al. (2009) have likely larger photometric errors 



than in our case, because their photometry comes from a short 
(2s) and a long (15s) FORS2-VLT image in V and / bands. 
Their larger photometric errors make difficult the identification 
of fSGB stars especially in the lower part of the SGB where the 
faint and the bright SGB have a smaller color difference. 

The CMD in the right panel of Fig. [3] suggests that the RGB 
of NGC 1851 has some spread in the B - I color. This spread 
could be related to the presence of the two RGB branches 
observed by Calamida et al. (2007) in the Stromgren (oti, 
u - y) CMD which is possibly associated with the presence 
of two groups of stars with different CN abundances (Yong 
& Grundahl 2008). To investigate how the upper and lower 
SGB populations may vary with radius, we divided the clus- 
ter into seven concentric annuli. The inner three are covered by 
the ACS data set out to 2.5 arcmin, and the outer four are cov- 
ered by the ground-based data set, going out to 12.5 arcmin. 
In Fig.|5]we show the CMD of the region around the SGB for 
stars in each annulus. Both the fainter and the brighter SGBs 
are clearly visible for radial distances smaller than ~8 arcmin. 
Beyond this radius, the number of stars is simply too small to 
distinguish any feature of the SGB. The WFPC2 observations 
allow an additional check on the distribution of stars along the 
SGB in the interval between 0.8 and 3.5 arcmin (see the central 
panel in Fig.|5]). 

4. Radial distribution of the population ratio 

4.1. The SGB subpopulations 

To determine the fraction of fSGB and bSGB stars, we adopted 
a procedure similar to that used by Milone at al. (2009). Fig.|6] 
illustrates this four-step procedure for the ACS/WFC sample. 

We selected by hand two points on the fSGB (P\j,P2j) and 
two points on the bSGB (P\^b,P2,b) with the aim of delimiting 
the SGB region where the split is most evident. These points 
define the two lines in panel (a), and only stars contained in the 
region between these lines were used in the following analysis. 

In panel (b) we have transformed the CMD linearly (the 
transformation equation is given in Appendix) into a reference 
frame where: the origin corresponds to P] j,; Pi j is mapped into 
(1,0), and the coordinates of Pjj, and Pjj are (0,1) and (1,1) 
respectively. For convenience, in the following, we indicate as 
'abscissa' and 'ordinate' the abscissa and the ordinate of this 
reference frame. The dashed green line is the fiducial of the 
bSGB. We drew it by marking several points on the bSGB, and 
interpolating a line through them by means of a spline fit. The 
black lines of panel (a) correspond to the loci with 'abscissa' 
of zero and one and to the loci with 'ordinate' of zero and one. 

In panel (c) we have calculated the difference between the 
'abscissa' of each star and the 'abscissa' of the fiducial Une (A 
'abscissa'). 

The histograms in panel (d) are the distributions in A 'ab- 
scissa' for stars in four A 'ordinate' intervals. These distribu- 
tions have been modeled as the sum of two partially overlap- 
ping Gaussian functions. To reduce the influence of outliers 
(such as stars with poor photometry, field stars and binaries) 
we did a preliminary fit of the Gaussian distribution using all 
available stars. Then we rejected all the stars more than two cr/, 
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F450W '"FS14»' 



Fig. 2. CMD of NGC 1851 from ACS (left) and WFPC2 (right) data. The inset show a zoom around the SGB region. 

NGC1851, ground-based data, 2.5<R<il.7 arcmin NGC1851, ground-based data, 2.5<R<11.7 arcmin 
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Fig. 3. y vs. V-I (left) and V ws. B - I (right) CMD of NGC 1851 from ground-based data. The inset show a zoom around the 
SGB region. Only stars with radial distance greater than 2.5 arcmin are plotted. 



to the left of the bSGB and more than two cr^ to the right of the 
fSGB and repeated the fit (the cr's are those of the best fitting 
Gaussian in each A'ordinate' bin fitted to the fSGB and bSGB 
respectively). In panel (c) the continuous vertical lines indicate 
the centers of the best-fitting Gaussians in each A'ordinate' in- 
terval. The red dashed line is located two crb on the left side of 
the bSGB, and the blue dashed line runs two crj on the right 
side ofthefSGB. 

We repeated the procedure for the ground-based sample, as 
illustrated in Figure [T] In this case we have reduced the num- 
ber of bins in 'magnitude' to just two intervals, because of the 
smaller number of stars. The same procedure has also been ap- 



plied to the WFPC2 data. In this case, we use the entire data 
set, without splitting it into bins, as illustrated in Fig. [8] 

It is important to notice that each of the points f i ^, Pij, 
P2,b, and f 2,/ — which we have arbitrarily defined with the sole 
purpose of isolating a group of stars representative of each of 
the two SGBs — corresponds to a different mass ( Mp\h, Mp\ /, 
Mp2h, and Mpif)- 

To obtain a more accurate measure of the fraction of stars 
in each of the two populations (hereafter: /bscB; /fSGs) we 
have to compensate for the fact that the two stellar groups 
that define the two SGBs cover two different mass intervals 
(Mp2 f - Mpif Mp2b- Mp\b), due to the different evolution- 
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Fig. 4. y vs. B - / CMD of NGC 1851 from ground-based data for stars with radial distances from the cluster center larger than 
2.5 arcmin (left) and CMDs the four quadrants right. 



ary lifetimes. Consequently, the correction we have to apply 
will be somewhat dependent on the choice of the mass func- 
tion. 

To this end, we can calculate the fraction of stars in each 
branch as; 



/bSGB - 
/fSGB - 



f Ni,INf 

Af 

NbINf 



where At, and Af are the areas of the Gaussians that best fit the 

bSGB and the fSGB, and Nf(b) = ^p'''"'' (f>(M)dM. 

As for the dependence on the adopted mass function, we 
ran the following test. Because of the effect of mass segrega- 
tion, we assumed a heavy-mass-dominated mass function for 
the central regions (a = -1.0), and a steep (a - 3.35) mass 
function for the external regions. Even with these extreme as- 
sumptions, we found that the mass function effect can change 
the relative fSGB/bSGB population ratio by a negligible 4%, 
in the small mass interval covered by our SGB stellar groups. 
Therefore, for simplicity, we adopted a Salpeter (1955) IMF for 
4>(M). 

Since the inferred population ratio depends on the as- 
sumed evolutionary lifetimes, our result is necessarily some- 
what model-dependent. In fact, the values could even be sce- 
nario dependent as well. Cassisi et al. (2008) propose that the 
two SGBs in NGC 1851 could in principle be explained by one 
of the following scenarios: 



- A) the bSGB and the fSGB are both populated by two nor- 
mal a-enhanced stellar populations with ages of 10 and 1 1 
Gyrs, respectively; 

- B) the fSGB belongs to an 1 1 Gyr old, normal a-enhanced 
stellar population, while the bSGB belongs to a population 
2 Gyr younger with high C+N+O abundances; 

- C) the bSGB and the fSGB both coiTespond to a 10 Gyr 
old, normal a-enhanced stellar population, but the fSGB 
population has C+N-i-O abundances a factor of two higher 
than the bSGB population. 

Salaris et al. (2008) found that scenario C does the best 
job of reproducing both the spectroscopic and photometric ob- 
servations, while scenario B seems the least consistent with 
the data. Regardless, these three possibilities illustrate that in- 
ferring the actual frequency ratio of the populations depends 
upon the assumed astrophysical explanation of the difference 
between them, through the different implications for the evolu- 
tionary lifetimes between our empirically chosen fiducial evo- 
lutionary states. 

To better characterize the radial distribution of SGB stars 
we divided both the ACS and ground-based data into concen- 
tric annuli as indicated in Fig.[T] The radial intervals have been 
chosen so that the each of the four ACS bins have the same 
number of stars, and the three ground-based bins also have the 
same number of stars among themselves. In each of these an- 
nuli, we calculated the fraction of bSGB and fSGB stars using 
the procedure described above. Table|2]summarizes the results. 
The first two columns list the radial interval and the average 
radial distance for SGB stars within this radial interval (Rave)- 
The origin of the data and the number of SGB stars in each bin 
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Fig. 5. CMDs of the region around the SGB of NGC 1851 at different radial distances, from the center out to the tidal radius, 
and from the different data bases. The double SGB is visible out to where there are enough stars to see the SGB. 



are indicated in the third and forth columns, respectively, while 
the last three columns list the ratio of faint to bright SGB stars 
that we obtain for each of the three scenarios of Cassisi et al. 
(2008). 

Since the values of the subpopulation ratios depend on the 
mass interval corresponding to each SGB segment, the results 
differ slightly depending on whether we assume that the SGB 
split corresponds to scenario A, B, or C, as we use stellar 
masses from three different sets of isochrones. The main re- 
sult of the paper is summarized in Fig.|9] Although Scenarios 
A and B show a hint of a slight radial gradient, it is not statis- 



tically significant. Scenario C does not even show a hint of a 
gradient. Even when we maximize our statistics by putting the 
ACS data in one bin, the WFPC2 photometry in a second bin, 
and all the ground-based data in a third bin (see first three rows 
in Table 2), there is no significant variation in the fSGB:bSGB 
ratio with cluster radius. 

When comparing the three different scenarios, we note that 
the radial trends are slightly different. This has no reason to be, 
since different evolutionary times should not change the rela- 
tive trends of the ratios. The problem comes from the fact that 
we are using three different data-sets, and we can not avoid 
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NGC 1851 ACS/WFC data R<2-5 arcmin 
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-2 2 
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-I 1 
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Fig. 6. This figure illustrates the procedure adopted to measure the fraction of stars belonging to the bSGB and fSGB in 
NGC 1851. Panel (a) shows a zoom of the ACS/WFC CMD. The two lines delimit the portion of the CMD where the split 
is most evident. Only stars from this region are used to measure the population ratio. In Panel (b) we have transformed the 
reference frame of Panel (a). The green dashed line is the fiducial of the region around the bSGB. In Panel (c) we plotted stars 
between the two lines but after the subtraction of the fiducial line 'abscissa' The four right bottom panels show the A'abscissa' 
distribution for stars in four A'ordinate' bins. The sohd lines represent a bigaussian fit. For each bin, the dispersions of the best 
fitting Gaussian are indicated. 



small transformation errors from theoretical to observational 
plane. ACS, WFPC2, and ground-based data require slightly 
different corrections because we are working with three dififer- 
ent filter systems, and we have used the models to infer the rel- 
evant mass ranges. These uncertainties generate small system- 
atic errors in our estimate of the fSGB:bSGB ratio. However, 



these errors are well within our estimated error bars, and do not 
change the main result of a flat trend. 

We want to emphasize that the fractions of faint to bright 
SGB stars calculated in this paper and listed in Table 2 accu- 
rately account for the mass interval covered by each SGB seg- 
ment. Therefore they differ from the numbers quoted in Milone 
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Fig. 7. As in Fig.|6]for stars from the ground-based sample and with 2.5<R<8.0 arcmin. 
Table 2. Percentage of fSGB and bSGB stars at several radial distances. 
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Fig. 8. As in Fig.|6]for stars with 0.8<R<3.5 arcmin. In this case we used the WFPC2 photometry. 



et al. (2008) where we have simply determined the number of 
stars belonging to two segments of fSGB and bSGB in a given 
wipeoew — '«F8i4w color interval. 

4.2. The HB subpopulation 

As we mentioned above, NGC 1851 is a prototypical bimodal 
HB cluster. In Fig. [10] we show the HB region for several an- 
nular bins both from ACS/WFC (top panels) and ground-based 
(bottom panels) data. Red HB stars are marked with red sym- 
bols, while blue HB stars and RR Lyrae are indicated in blue. 



Table 3 lists the ratio between the blue HB stars plus RR 
Lyrae and red HB stars in three different radial annuli. In calcu- 
lating these numbers we were careful to account for the fact that 
the lifetime of star in the blue HB is on average 11 % greater 
than the lifetime in the red HB. This estimate has been obtained 
by comparing the core He-burning lifetime for the mean stel- 
lar mass populating the red HB with that of the mean stellar 
mass populating the blue HB in the HB synthetic models made 
by Salaris et al. (2008) for simulating the HB distribution in 
NGC1851. 

The small number of intervals into which we have divided 
the field of view is a consequence of our need for a statist!- 



10 



Milone et al.: NGC 1851 SGB's gradient 




I [ I I li ] I I I I I I I I I I I i"1 r I I I I li I I I I I I I I I I I I I 
1.5 2 2.5 1.5 2 2.5 

log (RADIUS) 

-• — ACS data " WFPC2 data gB data 



Fig. 9. Fraction of fainter over brighter SGB stars for the scenarios A, B, and C and fraction of blue HB over red HB. Circles, 
squares and triangles refer to the ACS/WFC, WFPC2 and WFl data sets. The dashed vertical lines mark the core and the half-mass 
radius. 



cally significant number of stars in each subsample. The bot- 
tom right panel of Fig. |9] shows the trend of the (blue HB + 
RR Lyrae):(red HB) ratio as a function of the radial distance 
from the cluster center. 

Note that for these bright objects we can now safely include 
the ground-based photometry within 2.5 arcmin from the clus- 
ter center, going as close as 0.8 arcmin, before crowding con- 
ditions become prohibitive even for these luminous stars. The 
overlap regions of the two data-sets, between 0.8 and 2.5 ar- 
cmin, provide us with an important cross-check on the solidity 
and consistency of the two independent sets of photometry. 

Fig. |9] shows that there is no statistically significant evi- 
dence of a radial gradient in either the fSGB:bSGB ratio or the 
bHB:rHB ratio. Moreover, the results shown in Fig. 9 provide 
further support to the suggestion of Milone et al. (2008) and 
Salaris et al. (2008) that the blue HB stars are the progeny of 
the fSGB population, whilst the bulk of the red HB component 
is related to the bSGB population. More accurate spectroscop- 
ical measurements are mandatory in order to fully assess this 
scenario. 
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Table 3. Percentage of fSGB and bSGB stars at several radial 
distances. 
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Fig. 10. Zoom of the CMD region around the HB at several 
radial distances from ACS (top) and WFI data (bottom). We 
indicated in red the probable red HB stars and RR Lyrae and in 
blue the probable blue HB stars. 
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5. Discussion 

D'Ercole et al. (2008) have shown that if a second generation 
of stars is formed by material coming from or polluted by a 
first generation, then we would expect these stars to be born 
in the core of a stellar cluster, where a cooUng flow collects 
the gas ejected by the earlier population. As the cluster evolves 
dynamically, the two populations mix and the ratio of second 
over first generation tends to a constant value in the inner part 
of the cluster. Until mixing is complete, the radial profile of this 
ratio is fiat in the inner part and decreases in the outer cluster 
regions. By studying the radial profiles of the difi'erent popula- 
tions, we might hope to see evidence of these initial gradients 
before dynamical relaxation washes them out. 

Before this paper and the work by Zoccali et al. (2009), 
ci) Centauri was the only cluster where the radial distribution 
of different stellar sub-populations had been analyzed. In cj 
Centauri the stellar population associated with the blue, more 
metal rich MS is more centrally concentrated than the red, more 
metal poor one, with the relative ratio of blue over red MS star 
counts being quite constant within the cluster core beyond ~ 
12 arcmin and increasing by a factor of two from 8 arcmin 
(Solhma et al. 2007) to about 5 arcmin from the center. Then 
the blue MS to red MS ratio remains constant in the cluster core 
(BelHni et al. 2009). 

In this paper, we have investigated the radial distribution of 
the two stellar populations associated with the double SGB of 
NGC 1851. By couphng HST mA ground-based data we fol- 
lowed the distribution of the two populations from the center 
out to the tidal radius, both on the SGB and the HB. Salaris 
et al. (2008) claimed that the ratio of the two SGBs is con- 
sistent with the idea that the the progeny of the fSGB si dis- 
tributed from the blue to the red HB, including the whole insta- 
bility strip, while the bright one should be confined to the red. 
According to the scenario of Salaris et al. (2008), the fraction 
of fSGB that evolves into the red HB corresponds only to the 
5% of the total number of stars. 

At variance with Zoccali et al. (2009), we have clearly de- 
tected both the brighter and the fainter SGB at all radial dis- 
tances, out to ~8 arcmin. At larger radii, the number of SGB 
stars is simply too small to detect any substructure in the SGB. 
We have determined the ratio of fainter to brighter SGB stars 
and, unlike the case with co Centauri, found that — within the 
error bars — the two stellar populations have the same radial 
distribution. This conclusion does not depend on which of the 
three scenarios proposed by Cassisi et al. (2008) we assume to 
explain the observed SGB dichotomy. The ratio of the blue HB 
stars to the red HB stars also shows no significant trend with 
cluster radius. 

Whatever the origin and formation process of the two stel- 
lar generations in NGC 1851, now the two groups seem to be 
well mixed within the cluster. Because of the short relaxation 
time of NGC 1851 (log tyh ~ 8.8) this result may be not to- 
tally unexpected, as the cluster stars must be at least partially 
mixed. Unlike NGC 1851, co Centauri has a very long relax- 
ation time (log fr/i ~ 10), and it must have retained some in- 
formation on the radial distribution of the gas from which its 
multiple stellar generations formed. However, only a detailed 
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dynamical model, following the suggestion of D'Ercole et al. 
(2008) and Decressin et al. (2007), will answer the question 
of whether the gradient seen in Omega Cen and the lack of a 
gradient seen in NGC 1851 imply different initial radial distri- 
butions for the various populations. 
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Appendix Upon request from the referee, we give here the 
transformation equations used to pass from panel (a) to panel 
(b) in Figs. 6, 7, and 8. For simplicity we indicate the color and 
the magnitude as C and M respectively, and the 'abscissa' and 
the 'ordinate' as abs and ord. 

abs = (CR - c2)/c3 (I) 
ord = yr/cl (2) 
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where: 
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CR = {C- Cpi,b) cose + (M - Mpi,6) sine (3) 

MR = -(C - Cpi,6) sine + (M - Mpib) cose (4) 

MRpi^h - MRp2,b 

= atan—— — (5) 

CRpi^b - CRp2 b 

MRp2^b - MRp2j CRp2,bMRp2f - CRp2fMRp2,b 

cl = -CR+ (6) 

CRp2,b — CRp2j CRp2,b — CRp2,f 

CRpx^b - CRp2,b , ordpij,CRp2,b - ordp2,bCRpi,b 

= ~5 1 — ^'"^ + 1 1 

ordpij, - orap2,b orap\^ - orap2,b 

CRpij - CRp2,f - c2pi,f + clp2j ordpij(CRp2j - c2p2j) - ordp2.f(CRpij - clpij) 

ci ordH (0) 

ordpxj - ordpij ordp\j - ordp2j 



